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Abstract 

This is a survey of recent studies of singularity formation in solutions of spheri- 
cally symmetric Yang-Mills equations in higher dimensions. The main attention is 
focused on five space dimensions because this case exhibits interesting similarities 
with Einstein's equations in the physical dimension, in particular the dynamics 
at the threshold of singularity formation shares many features (such as universal- 
ity, self-similarity, and scaling) with critical phenomena in gravitational collapse. 
The borderline case of four space dimensions is also analyzed and the formation 
of singularities is shown to be intimately tied to the existence of the instanton 
solution. 

1 Introduction 

One of the most interesting features of many nonlinear evolution equations is the spon- 
taneous onset of singularities in solutions starting from perfectly smooth initial data. 
Such a phenomenon, usually called "blowup", has been a subject of intensive studies 
in many fields ranging from fluid dynamics to general relativity. Whether or not the 
blowup can occur for a given nonlinear evolution equation is the central mathematical 
question which, from the physical point of view, has a direct bearing on our understand- 
ing of the limits of validity of the corresponding model. Unfortunately, this is often a 
difficult question. Two famous examples for which the answer is not known are the 
Navier-Stokes equation and the Einstein equations. Once the existence of blowup is 
established for a particular equation, many further questions come up, such as: When 
and where does the blowup occur? What is the character of blowup and is it universal? 
Can a solution be continued past the singularity? 

In this paper we consider these questions for the Yang-Mills (YM) equations in higher 
dimensions. In the physical 3 + 1 dimensions, where the YM equations are the basic 
equations of gauge theories describing the weak and strong interactions of elementary 
particles, it is known that no singularities can form. This was shown by Eardley and 
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Moncrief |I[ who proved that solutions starting from smooth initial data remain smooth 
for all future times. The motivation for studying the YM equations in unphysical D + l 
dimensions for D > 3 is twofold (and unrelated to the latest fashion of doing physics 
in extra dimensions). From the mathematical point of view, it is the obvious thing to 
ask how the property of global regularity depends on the dimension of the underlying 
spacetime and whether singularities can form in D + 1 dimensions for D > 3. However, 
there is also a less evident physical reason which is motivated by the hope that by 
understanding the problem of singularity formation for the YM equations one might 
get insight into the analogous, but much more difficult, problem in general relativity. 
From this viewpoint - in which the YM equations are considered as a toy model for 
the Einstein equations - it is essential that these two equations belong to the same 
criticality class. Let us recall that the criticality class is defined as the degree a in the 
homogeneous scaling of energy E —* X a E under dilations x — > x/X. The classification 
of equations into subcritical (a < 0), critical (a = 0), and supercritical (a > 0) is a 
basis of the heuristic meta-principle according to which subcritical equations are globally 
regular, while supercritical equations may develop singularities for some (large) initial 
data 0. For the YM equations we have olym — D — 4, while for the Einstein equations 
oiE = D — 2. Therefore, the YM equations in D = 5 have the same criticality, a = 1, as 
the Einstein equations in the physical dimension. Another way of seeing this is to note 
that in D = 5 the dimension of the YM coupling constant [e 2 ] = M _1 L D_4 (in c = 1 
units) is the same as the dimension of the physical Newton's constant [G] = M~ X L. 

For the reason just explained, the main body of this paper is focused on the lowest 
super-critical dimension D = 5. In Section 3 we show that in this case there exists a 
countable family of regular (by regularity we mean analyticity inside the future light 
cone) spherically symmetric self-similar solutions labelled by a nonnegative integer n 
(a nodal number). Next, using linear stability analysis we show in Section 4 that the 
number of unstable modes around a given solution is equal to its nodal number. The 
role of self-similar solutions in the dynamical evolution is studied in Section 5, where 
we show that: i) the n = solution determines a universal asymptotics of singularity 
formation for solutions starting from generic "large" initial data; ii) the n — 1 solution 
plays the role of a critical solution sitting at the threshold of singularity formation. The 
latter is in many respects similar to the critical behaviour at the threshold of black 
hole formation in gravitational collapse. In both cases the threshold of singularity (or 
black hole) formation can be identified with the co dimension-one stable manifold of a 
self-similar solution with exactly one unstable mode. These similarities are discussed in 
detail in Section 6. 

We consider also the Cauchy problem for the YM equations in D = 4. Despite 
intensive studies of this borderline case, the problem of global existence is open. In 
Section 7 we describe numerical simulations which, in combination with analytic results, 
strongly suggest that large-energy solutions do blow up. We show that the process 
of singularity formation is due to concentration of energy and proceeds via adiabatic 
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shrinking of the instanton solution. At the end, a recent attempt of determining the 
asymptotic rate of shrinking is sketched. 

We remark that there are close parallels between YM equations in D + 1 dimensions 
and wave maps in (D — 2) + 1 dimensions 0. Indeed, many of the phenomena described 
here are mirrored for the equivariant wave maps into spheres in three [§], [5[] and two |J 
spatial dimensions. 

Sections 5 and 7 of this survey are based on joint work with Z. Tabor [[/J. The 
material of Sections 3 and 4 is new. 



2 Setup 

We consider Yang-Mills (YM) fields in D + 1 dimensional Minkowski spacetime (in the 
following Latin and Greek indices take the values 1,2, ... ,D and 0, 1, 2, . . . , D respec- 
tively). The gauge potential A a is a one-form with values in the Lie algebra g of a 
compact Lie group G. In terms of the curvature F a p = d a Ap — dpA a + [A a , Ap] the 
action is 



S = — j Tr(F a pF af5 )d D x dt, 



where e is the gauge coupling constant. Hereafter we set e = 1. The YM equations 
derived from ([!]) are 

d a F a ? + [A a ,F a P] =0. (2) 

As written, this equation is underdetermined because of the gauge invariance 

A a -> U^AJJ + U- l d a U, (3) 

where U is an arbitrary function with values in G. In order to correctly formulate 
the Cauchy problem for equation (0), one must impose additional conditions which fix 
this gauge ambiguity. We shall not discuss this issue here because in the spherically 
symmetric ansatz, to which this paper is restricted, the gauge is fixed automatically. 

For simplicity, we take here G = SO(D) so the elements of so(D) can be considered as 
skew-symmetric Dx D matrices and the Lie bracket is the usual commutator. Assuming 
the spherically symmetric ansatz || 

4{x) = {^-8y) l - w r ^ r \ (4) 

the YM equations reduce to the scalar semilinear wave equation for the magnetic po- 
tential w(t, r) 

D — 2 

-w tt + A( D _ 2) wH -g— w(l - w 2 ) = 0, (5) 
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where Am_ 2 ) = <9 2 + ^-^-d r is the radial Laplacian in D — 2 dimensions. The central 
question for equation (§) is: can solutions starting from smooth initial data 

w(0,r) = f(r), w t (0,r)=g(r) (6) 

become singular in future? As mentioned above, in the physical D = 3 dimensions Eard- 
ley and Moncrief answered this question in the negative Q. However, simple heuristic 
arguments indicate that the property of global regularity enjoyed by the YM equations 
in D = 3 might break down in higher dimensions. In order to see why the global be- 
haviour of solutions is expected to depend critically on the dimension D, we recall two 
basic facts. The first fact is the conservation of (positive definite) energy 

oo 

E = Tr {Fl + fg) d D x = c(D) J (w] + w 2 r + ^(l- W 2 ) 2 ) r D ~ 3 dr, (7) 

o 

where the coefficient c(D) — (D — l)fo/(S' D_1 ) follows from the integration over the 
angles and taking the trace. The second fact is scale-invariance of the YM equations: 
if A a (x) is a solution of (0), so is A a (x) = X^A^x/X), or equivalently, if w(t,r) is 
a solution of (|5|), so is w(t,r) = w(t/X,r/X). Under this scaling the energy scales as 
E = X D ~ A E, hence the YM equations are subcritical for D < 3, critical for D = 4, and 
supercritical for D > 5. In the subcritical case, shrinking of solutions to arbitrarily small 
scales costs infinite amount of energy, so it is forbidden by energy conservation. In other 
words, transfer of energy to arbitrarily high frequencies is impossible and consequently 
the Cauchy problem should be well posed in the energy norm. This important fact 
was proved in D = 3 by Klainerman and Machedon M, who thereby strengthened the 
result of Eardley and Moncrief. In the supercritical case, shrinking of solutions might 
be energetically favourable and consequently singularities are anticipated. In fact, we 
shall show below that singularities do form already in the lowest supercritical dimension 
D = 5. In the critical dimension D = 4 the problem of singularity formation is more 
subtle because the scaling argument is inconclusive. 



3 Self-similar solutions in D = 5 

In order to set the stage for the discussion of singularity formation we first need to 
analyze in detail the structure of self-similar solutions of equation (|5|). As we shall see, 
these solutions play a key role in understanding the nature of blowup. By definition, 
self-similar solutions are invariant under dilations w(t,r) — > w(t/X,r/X), hence they 
have the form 

w(t,r) = W( V ), V = 7 F L ~, (8) 
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where a positive constant T, clearly allowed by the time translation invariance, is intro- 
duced for later convenience. Note that for a self-similar solution we have 

hence the solution becomes singular at the center when t — > T (there is no blowup in 
the first derivative because regularity demands that W'(0) = 0). Thus, each self-similar 
solution W(r]) provides an explicit example of a singularity developing in finite time 
from smooth initial data. 

Substituting the ansatz (§) into (|5|) one obtains the ordinary differential equation 

D-3 {D - 5)r)\ jxrl , D-2 jxr2 



W" + + ±- p W + ^- =rW(l - W 2 ) = 0. (10) 

\ rj 1 — r) z J r} z (l — rj z ) 

As explained in the introduction, because of the expected connections with the Einstein 
equations, we are mainly interested in the lowest super-critical dimension D = 5. In 
this case equation ( TO ) reduces to 

w" + -w'+ 9/i 3 9N W (1 - W 2 ) = 0. (11) 
?7 r) z (l — r) z ) 

Although the similarity coordinate r\ is natural in the discussion of singularity formation, 
it has a disadvantage of not covering the region t > T, in particular it does not extend 
to the future light cone of the point (T, 0). For this reason we define a new coordinate 
x = 1/rj which covers the whole spacetime: the past and the future light cones are 
located at x = 1 and x — — 1, respectively; while the center r = corresponds to x = oo 
(for t < T) and x = — oo (for t > T). In terms of x equation ( pTTj) becomes 

(x 2 - l)W" + 3W(1 - W 2 ) = 0. (12) 

We first consider this equation inside the past light cone, that is for 1 < x < oo and 
impose the boundary conditions 

W(l) = and W(oo) = ±1, (13) 

which follow from the demand of smoothness at the endpoints. As we shall see below, 
once a solution to this boundary value problem is constructed, its extension beyond the 
past light cone can be easily done. 

To show that equation (|12"D admits solutions satisfying (|13D we shall employ a shoot- 
ing technique. The main idea of this method is to replace the boundary value problem 
by the initial value problem with initial data imposed at one of the endpoints and then 
adjusting these data so that the solution hits the desired boundary value at the second 
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endpoint. In the case at hand we shall shoot from x = 1 towards infinity. Substituting a 
formal power series expansion about x = 1 into (0) one finds the asymptotic behaviour 

W{x) = a(x - 1) - ^{x -If + O ((x - l) 3 ) , (14) 

where a is a free parameter determining uniquely the whole series. In the following a 
solution of equation (|1^) starting at x = 1 with the asymptotic behaviour ([TJ|) will be 



called an a-orbit. Without loss of generality we may assume that a > 0. We claim that 
there is a countable set of values {a n } for which the a„-orbits exist for all x > 1 and 
have the desired asymptotics at infinity (such orbits will be called connecting). The 
proof consists of several steps. 

Step 1 (Local existence). First, we need to show that a-orbits do in fact exist, that is, the 
series ( |i~4|) has a nonzero radius of convergence. Since the point x — 1 is singular, this 
fact does not follow from standard theorems. Fortunately, in [TIJ Breitelohner, Forgacs, 



and Maison have derived the following result concerning the behaviour of solutions of a 
system of ordinary differential equations near a singular point: 

Theorem [BFM]. Consider a system of first order differential equations for n + m func- 
tions u = (ui, ...,u n ) and v = (vi, ...,v m ) 

y^T = y^fi{y,u,v), y^- = -A»^ + y Vi gi{y,u,v), (15) 
ay ay 

where constants Aj > and integers > 1 and let C be an open subset of R n such 
that the functions f and g are analytic in the neighbourhood of y = 0, u = c, v = for 



all c E C . Then there exists an n-parameter family of solutions of the system such 
that 

Ui {y) = q + 0(y^), Vi {y) = 0(y"% (16) 

where Ui(y) and Vi(y) are defined for all c G C, \y\ < yo(c) and are analytic in y and c. 

We shall make use of this theorem to prove the local existence of a-orbits. In order 
to put equation ( p!2[ ) into the form ( jl6|) we define the variables 



and get 



W 

y = x-l, u(y) = W, v(y) = --W, (17) 

x — 1 



yv = -v + yf, yu =yf, f = — . (18) 

2 + y 



Since the function f(y,u,v) is analytic near y = for any u and v, according to the 
BFM theorem there exists a one-parameter family of local solutions such that 

u(y) = a + 0(y), v(y)=0(y), (19) 
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Transforming back to the original variables we obtain the behaviour of a-orbits. 

Step 2 (A priori global behaviour). It follows immediately from (|T2| ) that for x > 1 a 
solution cannot have a maximum (resp. minimum) for W > 1 (resp. < —1). Thus, 
once the solution leaves the strip \W\ < 1, it cannot reenter it (actually, such a solution 
becomes singular for a finite x). It is also clear that as long as \W\ < 1 the solution 
cannot go singular. To derive the asymptotics at infinity of a-orbits that stay in the 
strip \W\ < 1 we shall make use of the following functional 

Q{x) = l -{x 2 -l)W' 2 -\{l-W 2 ) 2 . (20) 

For solutions of equation ([12]) we have 

Q'( x ) = xW' 2 , (21) 

so Q(x) is monotone increasing. Now, we shall show that solutions satisfying \W\ < 1 
for all x > 1 tend to W = ±1 as x — > oo. To see this, first notice that for such solutions 
Q must be negative because if Q(xq) > for some xq > 1 then \W'\ is strictly positive 
for x > Xo so the solution must leave the strip \W\ < 1 in finite time. Since Q' > and 
Q < 0, it follows that Q has a nonpositive limit at infinity which in turn implies by (|21|) 
that lim^oo xW = and by (|20|) that lim^oo W exists. By L'Hopital's rule we have 
lim^oo x 2 W" = and using ( |T2"D again, we get that lim^oo W equals ±1 or 0. The 
latter is impossible because then Q(oo) = —3/4 is a global minimum contradicting the 
fact that Q increases. Thus, W(oo) = ±1. 

Step 3 (i) {Behaviour of a-orbits for small a). Rescaling w(x) = W(x)/a we get 

(x 2 -l)w" + 3w(l-aV) = 0, w(l) = 0, w'(l) = l. (22) 

As a —>■ 0, the solutions of this equation tend uniformly on compact intervals to the 
solution of the limiting equation 

(x 2 - l)w" + 3w = (23) 

with the same initial condition. This equation can be solved explicitly but for the 
purpose of the argument it suffices to notice that its solution, call it wl(x), is oscillating 
at infinity. Since W(x, a) ~ awi(x) up to an arbitrarily large x if a is sufficiently small, 
it follows that the number of zeros of the solution W(x,a) tends to infinity as a — > 0. 

(ii) (Behaviour of a-orbits for large a.) We rescale the variables, setting y = a(x — 1), 
w(y) = W(x) to get 

y(y + 2a)w" + 3w(l - w 2 ) = 0, w(0) = 0, w'(0) = 1. (24) 

As a — ► oo, the solutions of this equation tend uniformly on compact intervals to the 
solution w(y) = y of the limiting equation w" = 0. Thus, W(x, a) w a(x — 1) for large 
a and therefore the a-orbit crosses W = 1 for a finite x. 
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Step 4 (Shooting argument). We define the set 

A = {a | W(x, a) strictly increases up to some xq where W(xo, a) = 1}. (25) 

We know from Step 3 that the set Aq is nonempty (because the a-orbits with large a 
belong to it) and bounded below (because the a-orbits with small a do not belong to 
it). Thus ao = inf A exists. The solution W(x,a ) cannot cross the line W = 1 at a 
finite x because the same would be true for nearby solutions, violating the definition of 
a . Thus, < W(x, do) < 1 for all x and hence, by Step 2, lim^oo W(x, a ) = 1. This 

def 

completes the proof of existence of the nodeless self-similar solution W (x) = W(x, a ). 

Next, let us consider the solution with a = a® — e for small e > 0. By the definition 
of cto there must be a point x where this solution attains a positive local maximum 
W(xq) < 1 and since no minima are possible for < W < 1, it follows that there must 
be a point x\ > xq where W(xi, b) = 0. We shall show that Q(xi, a) > provided that e 
is sufficiently small. As argued above this would imply that the solution W(x, a) leaves 
the strip \W\ < 1 via W = — 1. From ([H]) we have 



xi W(x ) 

Q(xi) - Q{x ) = J xW' 2 dx = - J xW'dW. (26) 

xo 

In order to estimate the last integral note that for x > xq 

Q(x) - Q(xo) = \{x 2 - l)W' 2 - 2(1 - W 2 f + 2(1 - W 2 (x )) 2 > 0, (27) 



so x\W'\ > - w2 ) 2 ~ (1 - W^Oro)) 2 . Substituting this into (H) one gets 

W(x ) 

Q(x 1 )>-^(1-W 2 (x )) 2 + ^ J y/(l-W^-(l-W^x )) 2 dW. 



:28) 



The right hand side of this inequality is equal to a/2/3 for W(x ) = 1 so, by continuity, 
it remains strictly positive for W(xq) near 1. By taking a sufficiently small e we can 
have W(xq) arbitrarily close to 1, hence Q{x\) > which proves that a-orbits with 
a = do — e have exactly one zero. This means that the set A\ = {a \ W(x, a) increases 
up to some Xq where it attains a positive local maximum W(xq) < 1 and then decreases 
monotonically up to some x\ where W(x\) = —1} is nonempty. Let a\ = inf A\. By 
Step 3, a\ exists and is strictly positive. Using the same argument as above we conclude 
that the ai-orbit must stay in the region \W\ < 1 for all x, hence lim^oo W(x, a±) = —1. 

def 

This completes the proof of existence of the self-similar solution W\{x) = W(x, a\) with 
exactly one zero. 
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The subsequent connecting orbits are obtained by induction. We conclude that 
there exists a countable family of self-similar solutions W n (x) indexed by the integer 
n = 0,1, ... n which counts the number of zeros for x > 1. 

Remark. Since the sequence {a n } is decreasing and bounded below by zero, it has a 
nonnegative limit lim ri ^ 00 a n = a* > 0. If a* > 0, then the a*-orbit cannot leave the 
region \W\ < 1 for a finite x (because the set of such orbits is clearly open) hence it 
must be a connecting orbit with some finite number of zeros. But this contradicts the 
fact that the number of zeros of ct n -orbits increases with n. Hence, a* = 0. This implies 
that for any finite x, W n (x) goes to zero when n — > oo. 

We remark that the existence of the solution Wq was first shown by Cazenave, 
Shatah, and Tahvildar-Zadeh H via a variational method. 

The shooting technique is not only a powerful analytical tool; it is also an efficient 
numerical method of solving two-point boundary value problems. The numerical results 
produced by this method are shown in Table 1 and Figure 1. 
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Figure 1: The first four self-similar solutions W n {x). 

Surprisingly, it turned out that = 5/4 (with very good accuracy). This was a hint 
that the solution Wo has a simple closed form. Indeed, playing with the power series 
expansion (|14) we found that 



Wn(x) 



X' 



x 2 + 



3 ' 



(29) 
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n 





1 


2 


3 


4 


5 


a n 


1.25 


0.4813158 


0.1864517 


0.0722966 


0.02803703 


0.01087315 



Table 1: The shooting parameters of solutions W n for n < 5. 



Below we show an amusing calculation by Maple which helped us in finding this formula. 

> restart ; 

> with(DEtools) : 

> with(numapprox) : 

> ode:=(x~2-l)*diff (w(x) ,x$2)+3*w(x)*(l- w (x) ~2)=0; 

ode :=(x 2 - 1) w(x)) + 3 w(x) (1 - w(x) 2 ) = 

> ic:=w(l)=0,D(w)(l)=5/4; 

IC ;=w(l)=0, DH(l)^ 

> sol:=dsolve({ode,ic},w(x)) ; 

sol : = 

> sol_f ormal : =rhs (dsolve({ode , ic} , w(x) ,type=series) ) ; 

soLformal := 

^ U _ X ) - l l ( x - 1)2 + ^ U _ 1)3 + Hi. U _ 1)4 _ j^L / _ 1)5 + q// _ 1)6 ) 



> pade_sol :=pade(sol_f ormal, x=l, [2,2] ) ; 

|(x-l) 2 + |x-| 
padesol := -^j 



4 + ^+g(^-i) 2 



> sol :=simplif y(pade_sol) ; 



sol : = 5 



x 2 -l 
3 + 5x 2 



> subs (w(x)=sol , ode) ; 



(x 2 -l) (t&(5-^-— -^)) + 15 



x 2 -l 



,9 W (X 2 — l) 2 



fe2 V 3 + 5 x 2 ' 



3 + 5a; 2 



= 



> simplify (%) ; 



= 
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So far our analysis of self-similar solutions was restricted to the interior of the past 
light cone of the singularity. To show that the solutions W n represent genuine naked 
singularities, we need to extend them to the future light cone, that is to x = — 1. 
Fortunately, such an extension creates no problem because an a-orbit shot backwards 
from x = 1 cannot go singular before reaching x = — 1. This follows immediately from 
([12]) by observing that, in the interval —1 < x < 1, W(x) is concave down (resp. up) 
for W > 1 (resp. W < —1), hence W(x) remains bounded as x — > — 1 + . Moreover, the 
function Q(x) is negative and decreasing near x = —1, thus linx J .^_ 1 + Q(x) exists which 
implies in turn that c = linx r ^_ 1 + W(x) exists. Having that, the standard asymptotic 
analysis gives the following leading order behaviour for x — > — 1 + 

W{x) ~ c + -c(l - c 2 ){x + 1) ln(x + 1). (30) 

The singular logarithmic term in (|30|) can be eliminated by fine-tuning the shooting 
parameter a, however this is not expected to happen for the solutions W n (x) because 
in their construction the freedom of adjusting a was already used to tune away the 
singular behaviour for x > 1. We conclude that the self-similar solutions W n are C° 
at the future light cone and are analytic everywhere below it. The only (somewhat 
surprising) exception is the solution Wq which is analytic in the entire spacetime. 



4 Linear stability of self- similar solutions 

In this section we study the linear stability of self-similar solutions W n . This analysis 
is essential in determining the role of self-similar solutions in dynamics. We restrict 
attention to the interior of the past light cone of the point (T, 0) and define the new 
time coordinate s = — In \J (T — t) 2 — r 2 . Note that s — > oo when t — > T, and the lines 
of constant s are orthogonal to the rays of constant x. In terms of s and x, equation 
d^) becomes (for D — 5) 

e 2s 

-^—(e- 2s w s ) s + {x 2 - IK, + 3w(l - w 2 ) = 0. (31) 
x z — 1 

Of course, this equation reduces to (0) if w does not depend on s. In order to determine 
the stability of self-similar solutions W n we seek solutions of ([H]) in the form w(s,x) = 
W n (x) + v(s, x). Neglecting the 0(v 2 ) terms we obtain the linear evolution equation for 
the perturbation v(s,x) 

e 2s 

-^—{e~ 2s v s ) s + (x 2 - l)v xx + 3(1 - 3W 2 )v = 0. (32) 
x z — 1 

Substituting v(s,x) = e ( - a+1 - )s v / ^ 2 — 1 u(x) into (|32|) we get the eigenvalue problem in 
the standard Sturm-Liouville form 

~T ( {x2 ~ ^P) ~ 3(1 " 3W " )U = ~T~~T M ' ( 33 ) 
ax \ ax J x A — 1 
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where A = —a 2 . Using the variable p = §l n (fzj) ranging from zero to infinity we 
transform (|33"D into the radial Schrodinger equation 



d 2 u 
dp L 



+ V n u = Xu, 



V„ 



3(1 - sw 2 : 

sinh 2 p 



(34) 



The potential V n (p) has a typical "quantum mechanical" shape (see Figure 2) with the 
asymptotics 



V n (p) 



Q/p 2 for p — > 0, 
•12exp(— 2p) for p — > oo. 



(35) 



Note that the potential can be expressed in the form V n (p) = 1(1 + l)/p 2 + V^ eg (p) with 
I = 2, where the regular part V^ eg (p) is everywhere negative and V^ eg (0) — > — oo as 
n — > oo. 



3 




Figure 2: The potential for the perturbations around the self-similar solution W\. The single 
bound state with energy A = —16 is indicated. 



Both endpoints p = and p = oo are of the limit-point type, that is, exactly one 
solution near each point is square-integrable (admissible). Near p = the admissible 
solutions behave as u(p) ~ p 3 . For p — > oo and A < the admissible solutions behave 
as u(p) ~ e~ ap (recall that a = \J — A). All A > belong to the continuous spectrum. 

Let (resp. a£) denote the kth eigenfunction (resp. eigenvalue) about the solution 
W n . The numerically generated spectra are shown in Table 2. 
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n 


a 




a 2 


a 3 


a 4 
















1 





4 








2 





4 


27.407 






3 





4 


27.379 


182.49 




4 





4 


27.374 


182.18 


1214.5 


oo 





4 


27.37319 


182.1202 


1210.917 



Table 2: The eigenvalues of the perturbations about the first five solutions W n obtained 
numerically. The pseudo-eigenvalue a = is also included. The last row corresponding 
to n = oo was obtained by solving numerically the transcendental equation (ET3I). 



We point out that although a = is not a genuine eigenvalue, it is distinguished from 
the strictly positive part of the continuous spectrum by the fact that the corresponding 
non-square-integrable pseudo-eigenfunction, called the zero mode, is subdominant at 
infinity. The existence of the zero mode is due to the time translation symmetry, or in 
other words, the freedom of shifting the blowup time T in (§). To see this, consider the 
self-similar solution with a shifted blowup time W n ((T' — t)/r), where T' = T + e. In 
terms of the original similarity variables s = — In \/(T — t) 2 — r 2 and x = (T — t)/r, we 
have 



V -t 



W n (x + eeVx 2 - 1) = W n (x) + eeVx 2 - 1 W' n {x) + 0(e 2 



hence the perturbation generated by shifting the blowup time corresponds to a 
has the form 



(36) 



and 



Ur 



Vx 2 



1 WJx) 



sinh 2 pC(p)- 



(37) 
= 0, 



An alternative way of deriving this result is to take (sinh 2 p W')' + 3iy(l — W 2 ) 
which is (0) reexpressed in terms of p, differentiate it and compare with fl34l) . 

Since by construction the solution W n (p) has n extrema, it follows from (|37| ) that the 
zero mode Uq(p) has n nodes. This implies, by the standard result from Sturm-Liouville 
theory, that the potential V n has exactly n negative eigenvalues, in agreement with the 
numerical results shown in Table 2. We conclude that the self-similar solution W n has 
exactly n unstable modes (apart from the unphysical zero mode). In particular, the 
fundamental solution Wo is linearly stable, which makes it a candidate for the attractor. 

The rest of this section is a digression concerning a striking regularity which an 
acute reader might have already noticed in Table 2. Namely, the third column of 
Table 2 indicates that for each n > the first eigenvalue below the continuous spectrum 
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A™ = —(a™) 2 is equal to —16 (with the numerical accuracy of ten decimal places)! This 
puzzling numerical fact is calling for an explanation. Clearly, it has something to do 
with the particular form of the nonlinearity since, for instance, the analogous problem 
for self-similar wave maps from 3 + 1 dimensional Minkowski spacetime into the 3-sphere 
does not have this property Q . We suspect that the problem has some hidden symmetry, 
yet we cannot exclude a possibility that the numerics is misleading and the eigenvalues 
A" are not precisely equal but their splitting is beyond the numerical resolution. Some 
insight into this puzzle can be gained by analyzing the limiting case n — > oo. Recall 
that W n (p) tends to zero for any p > as n — > oo, hence the sequence of potentials V n 
has the following nonuniform limit 

lim V n (p) = = — ir- (38) 

n-»oo sinn p 

For the limiting potential the Schrodinger equation (|4]) can be solved exactly. The 
solution that is admissible at infinity (which as before is the limit-point) is given by the 
associated Legendre function of the first kind 



u(p) = P°(cothp), v = -1/2 + i^L. (39) 

Here v is one of the roots of v(y + 1) = —3; the second root gives the same solution 
because P-x/2+ip( x ) with real (3 is real and P^i/ 2+ ip( x ) = P-i/2-ip( x )- 
As p — > 0, the solution (|39|) behaves as 



v 11 

u(p) oc pa sin I - — In p + 5(a) 1 , (40) 

so it is always admissible, independently of a. This means that p = is the limit- 
circle point and therefore in order to have a well-defined self-adjoint problem we need 
to impose an additional boundary condition. In the language of spectral theory such 
a condition is called a self-adjoint extension. The continuous part of the spectrum is 
the same for all self-adjoint extensions but the eigenvalues do depend on the choice. In 
our case the self-adjoint extension amounts to fixing 5(a) - the phase of oscillations of 
the eigenfunctions for p — > 0. The natural choice is to require that the eigenfunctions 
oscillate with the same phase as the zero mode, that is 5(a) = 5(0), or equivalently 

lim (p°(cothpW(p) - P°'(cothpWp)) =0. (41) 

Note that under this condition the following diagram commutes 

v n * K } 



Ko K 00 } 
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Substituting (|39| ) into ( [41] ) and using the asymptotic expansion (for real (3) 

P« 1/2+ ^cothp) ~ ;^ +i/3 + C ' C - fOT ( 42 ) 

V27T r(l/2 + ip — a) 

we obtain the quantization condition for the eigenvalues 




iZ r )T{ l +i ^r +ak) \ = k7r > keN - (43) 

This transcendental equation has infinitely many roots which for k > 2 can be obtained 
only numerically (see the last row in Table 2). However, for k = 1 the exact solution is 
ax = 4 because (accidentally?) T(l/2 + + 4) = -45 T(l/2 + zVTT/2), as can be 

readily verified using four times the identity T{z + 1) = zT(z). Thus, we showed that 
the least negative eigenvalue of the limiting potential is equal to —16. Although this 
analysis does not resolve the original puzzle why all a" are equal to 4, it shows at least 
that 4 is the accumulation point of this sequence. 

The asymptotic distribution of eigenvalues for k — > oo can be derived from ( fl3| ) by 
using the formula for the asymptotic behaviour of the gamma function for large z 



r> + z) ~ V 7 ^ e {v+z - 1,2)Xnz - z for \z\ -> oo, (44) 

which yields 

-j — In a for a — >• oo. (45) 




Applying this to (|43|) one gets 

Qjfc+i 



2tt 



This formula was useful is providing an initial guess in the the numerical root finding 
procedure for equation (031) for large k. 



5 Singularities in D = 5 

Having learned about self-similar solutions, we are now prepared to understand the 
results of numerical studies, first reported in 0, of the Cauchy problem for the YM 
equation in five space dimensions 

2 3 

wtt = w rr + ~w r + — w(l - w 2 ). (47) 
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The main goal of these studies was to determine the asymptotics of blowup. Our 
numerical simulations were based on finite difference methods combined with adaptive 
mesh refinement. The latter were instrumental in resolving the structure of singularities 
developing on vanishingly small scales. We stress that a priori analytical insight into 
the problem, in particular the knowledge of self-similar solutions was very helpful in 
interpreting the numerical results. 

We solved equation fl4"7|) for a variety of initial conditions interpolating between 
small and large data. A typical example of such initial data is a Gaussian (ingoing or 
time-symmetric) of the form 

w(0,r) = 1 - Ar 2 exp [-cx{r - R) 2 ] , (48) 

with adjustable amplitude A and fixed parameters a and R. The global behaviour of 
solutions is qualitatively the same for all families of initial data and depends critically on 
the amplitude A (or any other parameter which controls the "strength" of initial data). 
For small amplitudes the solutions disperse, that is the energy is radiated away to 
infinity and in any compact region the solution approaches the vacuum solution w — 1. 
This is in agreement with general theorems on global existence for small initial data 0. 
Heuristically, this follows from the fact the for a small amplitude the nonlinearity is 
dominated by the dispersive effect of the linear wave operator. For large amplitudes we 
observe the development of two clearly separated regions: an outer region where the 
evolution is very slow and a rapidly evolving inner region where the solution attains a 
kink-like shape which shrinks in a self-similar manner to zero size in a finite time T. 
The kink is, of course, nothing else but the self-similar solution Wo (7^7)!] (see Figure 3). 

We summarize these findings in the following conjecture: 

Conjecture 1 (On blowup in D — 5). Solutions of equation ^% corresponding to 
sufficiently large initial data do blow up in finite time in the sense that w rr (t, 0) diverges 
as t y T for some T > 0. The universal asymptotic profile of blowup is given by the 
stable self-similar solution: 

lim«j(t, (T-t)r) = W (r). (49) 

t/T 

We think that the basic mechanism which is responsible for the observed asymptotic 
self-similarity of blowup can be viewed as the convergence to the lowest "energy" con- 



figuration. To see this, let us rewrite ( f47|) in terms of the similarity variable 77 and the 
slow time r = — ln(T — t) to get 



w TT +w T + 2r]w VT = (l - rf) (w m + -w v ) +^w(l - w 2 ). 

\ V J V 



(50) 



1 Throughout this section we use the similarity variable r\ = ijv—^ (rather than x) and abuse the 
notation by writing W n {r}) to denote W n (rj) — W n {x). 
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Figure 3: The upper plot shows the late time evolution of time symmetric initial data of the 
form (f4§|) with a = 10, R = 2, and A = 0.2. As the blowup progresses, the inner solution 
gradually attains the form of the stable self-similar solution Wo(r / (T — t)). The outer solution 
appears frozen on this timescale. In the lower plot the rescaled solutions w(t, (T — t))r are 
shown to collapse to the profile Wo(r) (solid line). 
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In this way the problem of blowup was converted into the problem of asymptotic be- 
haviour of solutions for r — > oo. The natural "energy" functional associated with this 
problem is 

l 

KW = J (^X + f*r^) d * ( 51 ) 



K(w) has a minimum at the self-similar solution Wo and saddle points with n unstable 
directions at solutions W n with n > 0. Since the wave equation (|50|) contains a damping 
term reflecting an outward flux of energy through the past light cone of the singularity, 
we suspect (but cannot prove) that K(w) decreases with time. If so, it is natural to 
expect that solutions will tend asymptotically to the minimum of K{w). 

We already know that solutions with small data disperse and solutions with large 
data blow up. The question is what happens in between. Using bisection, we found 
that along each interpolating family of initial data there is a threshold value of the 
parameter, say the amplitude A*, below which the solutions disperse and above which 
a singularity is formed. The evolution of initial data near the threshold was found to go 
through a transient phase which is universal, i.e. the same for all families. This inter- 
mediate attractor was identified as the self-similar solution W\. Having gone through 
this transient phase, at the end the solutions leave the intermediate attractor towards 
dispersal or blowup. This behaviour is shown in Figure 4 for the time-symmetric initial 
data of the form (f4"5|). 

The universality of the dynamics at the threshold of singularity formation can be 
understood heuristically as follows^. As we showed above, the self-similar solution W\ 
has exactly one unstable mode - in other words the stable manifold of this solution 
has codimension one and therefore generic one-parameter families of initial data do 
intersect it. The points of intersection correspond to critical initial data that converge 
asymptotically to W\. The marginally critical data, by continuity, initially remain 
close to the stable manifold and approach W\ for intermediate times but eventually are 
repelled from its vicinity along the one-dimensional unstable manifold (see Figure 5). 

According to this picture the universality of the nearly critical dynamics follows im- 
mediately from the fact that the same unstable mode dominates the evolution of all 
solutions. More precisely, the evolution of marginally critical solutions in the interme- 
diate asymptotics can be approximated as 

w(t, r) = Wx(rj) + c(A)(T - t)~ ai ~ l V(r)) + radiation, (52) 

where V\ is the single unstable mode with the eigenvalue a.\ = 4. The small constant 
c(A), which is the only vestige of the initial data, quantifies an admixture of the unstable 

2 This heuristic picture of the dynamics near the threshold, borrowed from dynamical systems theory, 
has been first given in the context of Einstein's equation - see Section 6 and fl^| . 
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Figure 4: The dynamics of time-symmetric initial data of the form ( [48|) with amplitudes that 
are fine-tuned to the threshold of singularity formation. The rescaled solution w(t, (T — t)r) 
is plotted against ln(r) for a sequence of intermediate times. Shown (solid and dashed lines) 
is the pair of solutions starting with marginally critical amplitudes A = A* ± e, where A* = 
0.144296087005405. Since e = 10 -15 , the two solutions are indistinguishable on the first 
seven frames. The convergence to the self-similar solution W\ (dotted line) is clearly seen 
in the intermediate asymptotics. The last two frames show the solutions departing from the 
intermediate attractor towards blowup and dispersal, respectively. 



mode - for precisely critical data c(A*) = 0. The time of departure from the intermediate 
attractor is determined by the time t* in which the unstable mode grows to a finite size, 
i.e., c(A)(T - f )~ Ql ~ 1 ~ 0(1). Using c(A) w c'(A*)(A - A*) and substituting a x = 4, 
we get T — t* ~ \A* — A] 1 ^. Various scaling laws can be derived from this. For example, 
consider solutions with marginally sub-threshold amplitudes A = A* — e. For such 
solutions the energy density 

, , wj w 2 r 3{l-w 2 ) 2 
e (^) = -T + 7f + 2r4 ( 53 ) 

initially grows at the center, attains a maximum at a certain time ~ t* and then drops 
to zero. Substituting (|52|) into (|53[) we get that e(t, 0) ~ (T — 1)~ 4 , and hence e(t*, 0) ~ 



-4/5 
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Figure 5: A schematic phase space picture of the dynamics at the threshold of singularity 
formation. 



6 Connection with critical phenomena in gravita- 
tional collapse 

The behaviour of solutions near the threshold of singularity formation described above 
shares many features with critical phenomena at the threshold of black hole formation 
in gravitational collapse. To explain these similarities, we now briefly recall the phe- 
nomenology and heuristics of the critical gravitational collapse. Consider a spherical 
shell of matter and let it collapse under its own weight. The dynamics of this process, 
modelled by Einstein's equations, can be understood intuitively in terms of the compe- 
tition between gravitational attraction and repulsive internal forces (due, for instance, 
to kinetic energy of matter or pressure). If the initial configuration is dilute, then the 
repulsive forces "win" and the collapsing matter will rebound or implode through the 
center, and eventually will disperse. On the other hand, if the density of matter is 
sufficiently large, some fraction of the initial mass will form a black hole. Critical gravi- 
tational collapse occurs when the attracting and repulsive forces governing the dynamics 
of this process are almost in balance, or in other words, the initial configuration is near 
the threshold of black hole formation. The systematic studies of critical gravitational 
collapse were launched in the early nineties by the seminal paper by Choptuik [TT 



in 



which he investigated numerically the collapse of a self-gravitating massless scalar field. 

Evolving initial data fine-tuned to the border between no-black-hole and black-hole 
spacetimes, Choptuik found the following unforseen phenomena near the threshold: 
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(i) universality: all initial data which are near the black hole threshold go through a 
universal transient period in their evolution during which they approach a certain inter- 
mediate attractor, before eventually dispersing or forming a black hole. This universal 
intermediate attractor is usually referred to as the critical solution. 

(ii) discrete self-similarity: the critical solution is discretely self-similar, that is it is 
invariant under dilations by a certain fixed factor A called the echoing period. 

(iii) black-hole mass scaling: for initial data that do form black holes, the masses of 
black holes satisfy the power law ~ e 7 where e is the distance to the threshold and 
7 is a universal (i.e., the same for all initial data) critical exponent. Thus, by fine tuning 
to the threshold one can make an arbitrarily tiny black hole. Put differently, there is no 
mass gap at the transition between black-hole and no-black-hole spacetimes. 

What Choptuik found for the scalar field, has been later observed in many other 
models of gravitational collapse, although the symmetry of the critical solution itself 
was found to depend on the model: in some cases the critical solution is self-similar 
(continuously or discretely), while in other cases the critical solution is static (or peri- 
odic). In the latter case black hole formation turns on with finite mass. These two kinds 
of critical behaviour are referred to as the the type II or type I criticality, respectively, to 
emphasize the formal analogy with second and first order phase transitions in statistical 
physics. We refer the interested reader to |L2| for an excellent review of the growing 
literature on critical gravitational collapse. 

The present understanding of critical behaviour in gravitational collapse is based on 
the same phase space picture as in Figure 5, that is, it is associated with the existence 
of a critical solution with exactly one unstable mode. This picture leads to some quan- 
titative predictions. In particular, in the case of type II critical collapse, an elementary 
dimensional analysis shows that the critical exponent 7 in the power law My, ~ e 7 is a 
reciprocal of the unstable eigenvalue of the critical solution. 

By now, the similarities between type II critical gravitational collapse and the dy- 
namics at the threshold of singularity formation in the 5 + 1 YM equations should be 
evident. This analogy, together with similar results for wave maps in 3 + 1 dimen- 
sions ||, 1 13 1, shows that the basic properties of critical collapse, such as universality, 
scaling, and self-similarity, first observed for Einstein's equations, actually have nothing 
to do with gravity and seem to be robust properties of supercritical nonlinear wave equa- 
tions. The obvious advantage of toy models, such as the one presented in this paper, is 
their simplicity which allowed to get a much better analytic grip on critical phenomena 
than in the case of Einstein's equations; in particular, it was possible to prove existence 
of the critical solution. The only characteristic property of type II critical collapse which 
so far has not found in simpler models (besides, of course, the absence of black holes 
which are replaced by singularities) is discrete self-similarity of the critical solution. It 
would be very interesting to design a toy model which exhibits discrete self-similarity at 
the threshold for singularity formation because this could give us insight into the origin 
of this mysterious symmetry. 
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7 Singularities in D = 4 



In this section we consider the Cauchy problem for the YM equation in four space 
dimensions 

w tt = w rr H — tiv H — 5-wfl — w 2 )- (54) 



We begin by recalling some facts concerning equation (j54|) which are be important in 
understanding the dynamics of singularity formation. First, we note that, in contrast 
to D — 5, there are no smooth self-similar solution in D = 4. This follows from the 
fact that in D dimensions the local solutions of equation flT0| ) near the past light cone 
behave as (1 — r/ 2 )^ - , hence they are not smooth if D is even (in particular, they are 
not different iable in D — 4). Although such singular self-similar solutions do exist, 
they cannot develop from smooth initial data and therefore they are not expected to 
participate in the dynamics. 

Second, D = 4 is the critical dimension in the sense that the energy (0) does not 
change under scaling. This means that, even though the model is scale invariant, a 
nontrivial finite energy static solution may exist^. In fact, such a static solution is well 
known 

w s(r) = \ = A- (55) 

1 + r z 

This is the instanton in the four- dimensional euclidean YM theory. Of course, by re- 
flection symmetry, — Wsif) is also the solution. Since the model is scale invariant, 
the solution Ws(r) generates an orbit of static solutions Wg(r) = Ws(r/X), where 
< A < oo. 

To analyze the linear stability of the instanton, we insert w(t,r) = Ws{r) + e lkt v(r) 
into ( |54"D and linearize. In this way we get the eigenvalue problem (the radial Schrodinger 
equation) 

d 2 1 d ,rt ^ ,2 ~\ t { n 2(1 -3WJ) 

+ V(r) v = k 2 v, V(r) = — ^. (56) 



dr 2 r dr J r 2 

This problem has a zero eigenvalue k 2 = which follows from scale invariance. The 
corresponding eigenfunction (so called zero mode) is determined by the perturbation 
generated by scaling 



Mr) = -^x W s« 



4r 

a=i (l+~ r 



rW s(r) = 77-^2- (57) 



3 Another way of seeing this is to notice that only in D = 4 the YM coupling constant e 2 provides 
the scale of energy. 
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Since the zero mode vo(r) has no nodes, it follows by the standard result from Sturm- 
Liouville theory that there are no negative eigenvalues, and eo ipso no unstable modes 
around Ws(r). Thus, the instanton is marginally stable. Note that the zero eigenvalue 
lies at the bottom of the continuous spectrum k 2 > 0, hence there is no spectral gap in 
the problem. 

After these preliminaries, we return to the discussion of the Cauchy problem for 
equation ([54]). For small energies the solutions disperse, in agreement with general 
theorems. For large energies, at first sight the global behaviour seems similar to the 
D = 5 before, near the center the solution attains the form of a kink which 

shrinks to zero size. However, this similarity is superficial because now the kink is not 
a self-similar solution (as no such solution exists). It turns out (see Figure 6) that the 
kink has the form of the scale- evolving instanton 



^^^ l ' (58) 

where a scaling factor X(t) is a positive function which tends to zero as t — > T. We 
summarize these findings in the following conjecture: 

Conjecture 2 (On blowup in D = 4). Solutions of equation with sufficiently 

large energy do blow up in finite time in the sense that w rr (t,0) diverges as t f T 
for some T > 0. The universal asymptotic profile of blowup is given by the instanton. 
More precisely, there exists a positive function X(t) \ for t /* T such that 

\imu(t,\(t)r) = W s (r). (59) 

t/T 

The key question which is left open in this conjecture is: what determines the evo- 
lution of the scaling factor A(i); in particular, what is the asymptotic behaviour of A(t) 
for t — > T? Numerical evidence shown in Figure 7 suggests that the rate of blowup goes 
asymptotically to zero, that is ( ' = d/dt) 

lim = - lim A = 0, (60) 

t^TT-t t^T ' v ' 

but it seems very hard to determine an exact asymptotics of \(t) from pure numerics^. 



Recently, an analytical approach to this problem has been suggested in ||15|| . Below 
we sketch the main idea of this approach. Let A4 = {Ws(r/\)\\ G M + } be a manifold 
of rescaled instantons (a one-dimensional center manifold). Assuming that a solution is 
in a neighbourhood of A4, we decompose it as 

w = W s (v)+v(t,r)), 7] = r/\{t). (61) 



This issue is also discussed by Linhart and Sadun in pA 
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Figure 6: The upper plot shows the formation of a singularity for large initial data of the form 
(48) with A = 0.5. The inner solution has the form of the scale-evolving instanton Ws{r / X(t)) 
with the scale factor A(i) going to zero slightly faster than linearly. In the lower plot the 
rescaled solutions are shown to collapse to the profile of the instanton Ws(r) (solid line). 
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Here v represents a small deviation of the solution from A4 and A is the collective 
coordinate on M.. To fix the splitting between these two parts we require that Pv = 0, 
where P is the projection on M.. Plugging fl6l]) into (|54D we get (' = d/drf) 



X 2 v - 2t]X\v +Lv + N(v) = XXt]Ws - X 2 (r] 2 W's + 2-qW' s ), (62) 
where L is the linear perturbation operator about the instanton 

_ _ d 2 Id 2(1 - 3W§) 



drf 1] drj i] 2 

and 



(63) 



N(v) = ^v 2 + -v\ (64) 

It is clear from (|62D that v = 0(X 2 ), hence for v to decay to zero as t — > T, the rate 
of blowup must go to zero as well. We stress this point to emphasize that the linear 
evolution of X(t), predicted for example by the geodesic approximation, is inconsistent 



with Conjecture 2. Next, by projecting equation ( p2|) on A4 and in the orthogonal 
direction, we get a coupled system consisting of a nonhomogeneous wave equation for 
v and an ordinary differential equation for A. Solving the first equation for v and 
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plugging the result into the second equation, we obtain in the lowest order the following 
modulation equation 

AA = | A 4 . (65) 
From this we get the leading order asymptotics for t — > T 

As shown in Figure 7 this result is in rough agreement with numerics. There are many 
possible sources of the apparent discrepancy. On the numerical side there are dis- 
cretization errors, an error in estimating the blowup time, or errors in computing A 
from the data. On the analytical side, there might be corrections to ( |66"D coming from 
the bounded region expansion and, more importantly, from the far field behaviour^. 
Finally, and in our opinion most likely reason of discrepancy is that the solution shown 
in Figure 7 has not yet reached the truly asymptotic regime and consequently the higher 
order corrections to formula (|66| ) are still significant. 

The issue of blowup rate is closely related to the problem of energy concentration in 
the singularity. To explain this, we define the kinetic and the potential energies at time 
t < T inside the past light cone of the singularity 

T-t T-t 

E K (t) = 6tt 2 J wfrdr, E P (t) = 6tt 2 J (w 2 r + ^ ~J ^ j rdr. (67) 
o o 



Substituting (58) into fl67]) we obtain 



T-t T-t 
Aft) Ap5 



E K (t) = 6vr 2 A 2 J W' s 2 rdr, E P (t) = 6tt 2 J (w' s 2 + (1 rdr. (68) 

o o 



Assuming flSTf), this implies that 



\imE K (t) = 0, \imE P (t) = 16tt 2 . (69) 

Thus, the energy equal to the energy of the instanton gets concentrated in the singularity. 
This means that in the process of blowup the excess energy must be radiated away from 
the inner region as the solution converges to the instanton. 

It is worth pointing out that the concentration of energy is a necessary condition for 
blowup in the critical dimension. To see this, suppose that the solution blows up at time 

5 The derivation of ( |65|) is not quite straightforward because of the presence of infrared divergencies 
which need to be regularized. 
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T and assume for contradiction that lim^T E(t) = 0. Then, by choosing a sufficiently 
small e > we can have E(T — e) arbitrarily small. This implies, by causality and global 
existence for small energy data, that the solution exists globally in time, contradicting 
the assumption. For radial equations of the type 

W tt = W rr + -W r H — , (70) 



a stronger result was proved by Struwe who basically showed that if the solution 



of ( |70"D blows up, then it must do so in the manner described in Conjecture 2. In this 
sense formation of singularities is intimately tied with the existence of a static solution. 

Finally, we address briefly the issue of the threshold of singularity formation. Using 
the technique described in Section 5, along each interpolating family of initial data we 
can determine the critical point separating blowup from dispersal. However, in contrast 
to the D = 5 case, we see no evidence for the existence of an intermediate attractor 
in the evolution of nearly critical data. This fact, together with a similar result for 
2 + 1 dimensional wave maps |6| suggests that in the critical dimension the transition 
between blowup and dispersal is not governed by any critical solution. We suspect that 
the evolution of precisely critical initial data still has the form (|58|) but the dynamics 
of the scaling factor is different than in fl66l) . This belief is based on the fact that in the 
evolution of marginally critical initial data we can clearly distinguish a transient phase 
during which X(t) drops very quickly and then, either reaches a minimum and starts 
growing (in the case of dispersal), or keeps decreasing with the "normal" rate (in the 
case of blowup). 



8 Conclusions 

There are two main lessons that we wanted to convey in this survey. The first lesson 
is that there are striking analogies between major evolution equations. In particular, 
the mechanism of blowup to a large extent is determined by the criticality class of the 
model. These analogies can be used to get insight into hard problems (such as singu- 
larity formation for Einstein's equations) by studying toy models which belong to the 
same criticality class. This approach is in the spirit of general philosophy expressed 
by David Hilbert in his famous lecture delivered before the International Congress of 
Mathematicians at Paris in 1900 [J17|| : "In dealing with mathematical problems, special- 
ization plays, as I believe, a still more important part than generalization. Perhaps in 
most cases where we seek in vain the answer to a question, the cause of the failure lies 
in the fact that problems simpler and easier than the one in hand have been either not 
at all or incompletely solved. All depends, then, on finding out these easier problems, 
and on solving them by means of devices as perfect as possible and of concepts capable 
of generalization. " 
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The second lesson is concerned with the interplay between numerical and analytical 
techniques. Accurate and reliable numerical simulation of singular behaviour is difficult 
and hard to assess. In order to keep track of a singularity developing on exceedingly 
small spatio-temporal scales, one needs sophisticated techniques such as adaptive mesh 
refinement. For these techniques the convergence and error analysis are lacking so 
extreme care is needed to make sure that the computed singularities are not numerical 
artifacts. For this reason, in order to feel confident about numerics it is important 
to have some analytical information, like existence of self-similar solutions. Without a 
theory, simulations alone do not provide ample evidence for the existence of a singularity. 
We believe that the interaction between numerical and analytical techniques, illustrated 
here by the studies of blowup, will become more and more important in future as we 
begin to attack more difficult problems. 
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